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I.  statement  of  Work 


One  of  the  Navy’s  most  important  weapon  systems  is  the  cruise  missle. 
Generally  a  missile  is  called  a  cruise  missile  if  its  speed  is  subsonic,  if  it 
uses  a  built-in  global  positioning  navigation  system  (GPS/INS),  and  if  its 
range  is  at  least  several  hundred  miles.  During  the  early  1990s,  cruise  mis¬ 
siles  launched  from  ships  and  submarines  were  used  with  great  success  in 
the  Persian  Gulf.  Targets  located  hundreds  of  miles  inland  were  success¬ 
fully  attacked  by  these  Navy  weapons  stationed  in  relatively  safe  positions 
offshore.  In  addition,  cruise  missiles  have  been  used  to  disable  SAM  sites, 
making  manned  aircraft  strikes  safer. 

A  mission  for  a  cruise  missile  is  defined  as  a  path  from  the  launch  area  to 
the  target.  The  job  of  the  mission  planner  is  to  determine  the  path  that 
the  missile  will  follow  and  program  the  missile  by  specifying  a  sequence 
of  waypoints  and  including  any  TERCOM  maps  that  are  available.  The 
missile  uses  its  dn-board  computational  facilities  and  its  GPS/INS  system 
to  guide  it  to  its  target. 

Missiles  are  never  launched  as  singles  and  are  launched  in  groups  to  form 
a  strike.  A  strike  plan  involves  developing  a  plan  for  many  missiles  all  of 
which  may  be  launched  during  a  short  period  of  time.  Frequently  two  or 
more  missiles  will  have  the  same  target  or  different  locations  within  the 
same  target.  For  example,  three  missiles  may  be  needed  to  hit  different 
areas  of  a  plant  to  incapacitate  the  plant. 

Cruise  missies  are  designed  to  be  low  flying  (terrain  following)  which  makes 
them  extremely  difficult  to  destroy  enroute  to  the  target.  The  mission  plan¬ 
ners  attempt  to  either  avoid  SAM  sites  or  fly  so  low  that  a  SAM  radar  system 
can  not  see  the  incoming  missile.  However,  cruise  missiles  are  vulnerable  to 
anti-aircraft  artillery  (triple  A)  if  a  hostile  ground  crew  is  given  sufficient 
warning.  If  the  triple  A  ground  crew  has  its  guns  aimed  at  the  correct  spot, 
a  missile  can  be  destroyed  by  a  single  round  striking  a  key  component. 


The  best  strategy  to  defeat  triple  A  is  to  vary  the  missions  of  the  cruise 
missiles.  If  the  first  missile  of  a  strike  passes  directly  over  a  triple  A  site, 
then  the  second  should  fly  a  few  hundred  yards  to  either  the  left  or  right. 

During  1996  our  research  team  develped  algorithms  designed  to  aid  a  mis¬ 
sion  planner  in  the  mission  development  process.  One  of  these  algorithms 
was  a  auto-router  based  on  geometric  heuristics  which  performed  well  on 
randomly-generated  test  cases.  In  this  research  we  present  (1)  an  enhance¬ 
ment  of  the  auto-router  algorithm  which  allows  greater  flexibility  in  the 
modelling  of  the  threat  regions  associated  with  SAM  sites,  and  (2)  a  new 
algorithm  designed  to  aid  a  strike  planner  in  the  development  of  multiple 
mission  plans  which  incorporate  variation  in  related  missions. 


II.  Path  Generation  with  Convex  Polygonal  Threats 


Using  current  technology,  missions  are  developed  using  a  two-step  process. 
In  the  first  step,  a  mission  planner  uses  a  two-dimensional  map  to  manually 
select  a  path  from  a  launch  site  to  the  target.  In  the  second  step  a  soft¬ 
ware  system  includes  the  vertical  dimension  and  develops  an  estimate  of  the 
probability  of  a  successful  mission.  This  process  may  be  repeated  several 
times  until  the  mission  planner  is  satisfied  with  the  mission  plan. 

In  1996  we  developed  an  auto-router  algorithm  to  aid  in  mission  planning 
which  represented  threat  regions  as  circles.  Within  this  framework,  the 
problem  of  developing  a  single  mission  for  a  single  cruise  missile  can  be 
viewed  as  a  computational  geometry  problem  in  which  we  seek  a  path  com¬ 
posed  of  line  segments  from  a  launch  site  to  a  target  site  which  skirts  se¬ 
lected  threat  regions.  The  basic  idea  on  which  the  algorithm  is  based  finds 
a  circumscribed  triangle  which  encloses  a  given  circle  (threat)  and  has  one 
vertex  at  a  current  missle  position  (waypoint).  A  missile  can  avoid  a  threat 
by  moving  along  the  sides  of  such  a  triangle.  In  empirical  tests  on  twenty 
randomly  generated  test  problems,  the  circumscribed  triangle  auto-router 
was  found  to  work  very  well.  This  basic  algorithm  was  enhanced  to  allow 
some  penetration  of  the  threat  regions  by  the  missle  path  while  satisfing  a 
probability  constraint  for  a  successful  traversal  based  on  proximity  to  the 
threats. 

According  to  analysts  at  the  Workshop  for  6.1  Research  on  Optimization 
held  in  Dahlgren,  Virginia  on  1-2  October  1996,  threats  can  come  in  all 
shapes  and  sizes.  When  routing  over  water,  missions  usually  avoid  islands 
(which  can  take  any  shape)  and  other  ships  which  look  more  like  ellipses 
than  circles.  To  accomodate  these  more  general  shapes  for  threat  regions 
we  have  enhanced  the  algorithm  by  allowing  the  mission  planner  to  model 
any  shape  defined  as  the  convex  hull  of  given  points. 


A  technical  report  giving  details  of  the  algorithm  resulting  from  our  de¬ 
velopment  work  appears  in  Appendix  A.  This  report  includes  a  qualitative 
evaluation  baaed  on  visual  displays  of  solution  paths  for  randomly  generated 
test  problems. 


III.  Multiple  Path  Generation  with  Separation 


The  auto-routers  we  have  developed  were  designed  only  to  aid  in  deter¬ 
mining  a  single  mission  for  a  single  missile.  A  strike  plan  involves  many 
missiles  and  many  targets  (25  missiles  is  not  unusual).  Also  the  probability 
of  a  safe  traversal  through  an  area  is  affected  by  multiple  missiles  being 
routed  through  an  area.  We  have  assumed  that  vulnerability  to  triple  A 
can  be  effectively  reduced  by  introducing  a  reasonable  amount  of  separa¬ 
tion  in  the  flight  paths  of  missiles  at  the  points  of  nearest  approach  to  threat 

regions. 

To  aid  a  strike  planner  we  have  developed  an  algorithm  which  allows  the 
specification  of  several  missle  location  and  target  pairs  together  with  the 
specification  of  the  number  of  multiple  missiles  to  be  launched  for  each  pair. 
For  each  pair,  a  mission  plan  is  developed  for  the  first  missle  by  applying  the 
auto-router.  Then  each  threat  region  approached  by  the  first  missle  mission 
plan  is  slightly  enlarged.  Then  a  mission  plan  is  developed  for  the  next 
missle  by  applying  the  auto-router  to  the  modified  configuration.  Again,  any 
threat  region  approach  by  the  next  missle  mission  plan  is  slightly  enlarged. 
Continuing  in  this  manner  assures  that  multiple  missies  will  have  some 
separation  when  approaching  threat  regions  along  similar  paths.  After  all 
mission  plans  have  been  generated  for  all  missies  associated  with  a  particular 
missle  location  and  target  pair,  all  threat  regions  are  reset  to  their  normal 
configuration. 

A  technical  report  giving  details  of  the  strike  planning  algorithm  resulting 
from  our  development  work  appears  in  Appendix  B.  This  report  includes 
visual  displays  of  solution  paths  for  twenty  randomly  generated  test  prob¬ 
lems.  Each  problem  has  five  missle  location  and  target  pairs  and  has  five 
missies  to  be  launched  per  pair. 
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Abstract 


This  manuscript  presents  mathematical  models  and  solution  algorithms  for 
the  problem  of  finding  a  path  composed  of  line  segments  from  a  given  missile 
location  to  a  given  target  location  in  the  presence  of  threats.  The  threats  are 
modelled  as  either  circles  or  convex  polygons.  The  paths  must  either  avoid 
the  threats  or  satisfy  some  probability  constraint  of  a  successful  traversal 
based  on  the  proximity  of  the  path  to  the  threats.  The  algorithms  are  all 
based  on  geometric  concepts  and  their  output  is  evaluated  qualitatively  based 
on  visual  displays  of  the  problem  and  solution  path. 
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1  INTRODUCTION 


During  recent  conflicts  in  the  Persian  Gulf,  the  U.S.  strategy  was  to  launch 
a  cruise  missile  strike  followed  by  manned  aircraft  strikes.  The  cruise  missile 
strike  would  involve  approximately  25  missiles  which  were  launched  from 
Navy  ships  and  submarines  stationed  in  relatively  safe  offshore  positions. 
Communication  systems  and  SAM  sites  were  frequently  targeted,  which  made 
any  following  manned  strikes  substantially  safer. 

A  cruise  missile  strike  consisting  of  25  missiles  requires  developing  a  mis¬ 
sion  for  each  missile.  A  mission  is  defined  as  a  path  composed  of  line  segments 
from  a  launch  site  to  a  target  site.  Using  current  technology,  missions  are 
constructed  using  a  two-step  process.  In  the  first  step,  a  mission  planner  uses 
a  two-dimensional  map  to  manually  select  a  path  (mission)  from  the  launch 
site  to  the  target.  A  software  system  takes  this  path,  adds  the  vertical  di¬ 
mension,  and  estimates  the  probability  of  success.  The  probability  of  success 
is  a  function  of  the  path  chosen  and  is  based  on  the  length  of  the  mission 
and  the  proximity  of  the  path  to  known  threats. 

These  weapons  are  designed  to  be  low  flying  (terrain-following)  which 
makes  them  extremely  difficult  to  destroy  en  route  to  the  target.  The  mis- 


planners  attempt  to  either  avoid  SAM  sites  or  to  fly  so  low  that  a  SAM 
radar  system  cannot  detect  the  incoming  missile.  However,  cruise  missiles 
are  vulnerable  to  anti-aircraft  artillery  (AAA)  if  the  ground  crew  is  given 
sufficient  warning.  If  the  AAA  ground  crew  has  their  guns  aimed  at  the 
correct  spot,  the  missile  can  be  destroyed  if  a  round  strikes  a  key  component. 
The  best  strategy  to  defeat  AAA  is  to  vary  the  missions  of  the  cruise  missiles. 

If  the  first  missile  of  the  strike  passes  directly  over  an  AAA  site,  then  the 
second  should  fly  a  few  hundred  yards  to  either  the  left  or  right. 

After  a=mission  plan  has  been  selected  for  a  missile,  the  mission  planner 
programs  the  missile  with  this  path  including  any  TERCOM  maps  that  will 
be  used.  A  mission  does  not  have  to  include  TERCOM  maps,  but  they  can 
be  used  as  part  of  the  navigation  system.  Some  missions  only  involve  GPS 

waypoints. 

Obviously,  a  graphical  system  which  allows  the  mission  planner  to  create 
mission  plans  with  a  point-and-click  system  will  be  an  improvement  over  a 
completely  manual  system.  However,  requiring  the  mission  planner  to  repeat 
this  task  25  times  to  create  a  strike  plan  will  be  tiresome.  Navy  analysts  indi¬ 
cate  that  a  better  system  would  provide  the  mission  planner  with  a  first-cut 
strike  plan  and  allow  the  planner  to  modify  some  of  the  missions  using  a 
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point-and-click  mechanism.  To  automatically  create  a  strike  plan  involving 
25  missiles  requires  25  mission  plans.  Preliminary  to  creating  such  a  strike 
plan,  we  must  be  able  to  generate  a  single  cruise  missile  mission.  The  objec¬ 
tive  of  this  investigation  is  to  develop  algorithms  to  automatically  generate 
a  single  cruise  missile  mission.  Based  on  the  results  of  this  manuscript,  we 
have  developed  an  algorithm  to  create  a  complete  mission  plan  (see  [6]). 


1.1  Description  of  the  Problem 

For  this  investigation,  the  threat  regions  are  modelled  as  convex  polygonal 
sets  or  circular  regions.  The  mission  plan  is  composed  of  a  set  of  connected 
line  segments,  and  adjacent  line  segments  have  a  restriction  on  the  maximum 
turn  angle.  The  mission  illustrated  in  Figure  1,  is  composed  of  three  line 
segments  iAB,BC,CD)  and  has  two  waypoints  {B,C).  The  turn  angles 
are  given  by  lEBC  and  IFCD.  Mathematically,  the  problem  of  finding  a 
feasible  mission  plan  may  be  stated  as  follows: 

Given  two  points  M  and  T  €  a  turn  angle  restriction  of  Z, 
and  K  convex  sets  Ri,...  ■, Rki  fiiid  a  set  of  line  segments 
[71,72],  [72,73],  . . .,  [7s,7s+i]  with  7i  =  M  and  7s+i  =  T 
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such  that  no  line  segment  intersects  the  interior  of  any  Ri,..., Rk 
and  no  turn  angle  exceeds  Z. 


Figure  1  about  here 


Let  r(M,  T,Z,Ri,...,  Rk)  or  simply  F  denote  the  set  of  feasible  missions 
for  a  given  missile.  Note  that  if  F  ^  then  there  are  an  infinite  number  of 
paths  in  F.  If  we  could  identify  the  objective  function  for  mission  planning, 
say  f(y)  for  ?/  €  F,  then  we  would  have  a  well-defined  optimization  problem. 
Unfortunately,  the  objective  function  for  this  problem  involves  some  unknown 
trade7ofF  between  the  length  of  the  path  from  M  to  T  and  the  number  of 
segments  needed  for  the  mission.  For  the  problem  illustrated  in  Figure  2, 
the  shortest  path  from  M  to  T  which  does  not  pass  through  the  threat  is 
the  line  segment  MA,  the  arc  AB,  and  the  line  segment  BT.  Approximating 
this  path  with  line  segments,  as  illustrated  in  Figure  2b  results  in  an  infinite 
number  of  such  segments.  Clearly,  mission  planners  are  not  interested  in 
plans  of  the  type  illustrated  in  Figure  2b.  For  the  problem  illustrated  in 
Figure  3,  a  two  segment  solution  is  given  by  MC,  CT.  However,  the  three 


segment  solution  MA.,AB,BT  is  shorter. 


Figures  2  and  3  about  here 


Unlike  other  optimization  problems,  the  mission  planning  problem  as 
we  describe  it  has  no  mathematical  objective  function.  The  objective  is 
to  find  short  feasible  missions  with  only  a  few  segments.  Given  a  feasible 
mission  plan,  we  believe  that  a  mission  planner  can  look  at  a  visual  display 

c 

of  the  mission  and  conclude  whether  or  not  this  is  a  good  plan.  Therefore, 
we  are  seeking  algorithms  which  find  feasible  missions  with  three  essential 
characteristics:  a  short  total  mission  distance,  a  small  number  of  segments, 
and  turn  angles  within  reasonable  limits. 

1.2  Survey  of  the  Literature 

Most  research  investigations  in  the  area  of  mission  planning  begin  by  placing 
a  grid  over  the  combat  theatre  and  applying  a  modification  of  Dijkstra  s 
algorithm  (see  [2])  to  obtain  a  safe  path  from  the  origin  to  the  destination 
(see  [1,  5, 16, 17,  27,  28]).  There  are  two  disadvantages  to  this  approach.  The 


grid  graph  can  become  very  large  and  consequently  the  algorithms  can  be 
very  time-consuming.  Also,  the  resulting  paths  may  involve  numerous  line 
segments,  some  of  which  may  have  90°  (excessive)  turn  angles.  This  approach 
uses  discrete  methods  for  a  problem  that  is  fundamentally  continuous.  We 
believe  that  continuous  methods  that  consider  the  geometry  of  the  problem 
are  preferable  to  discrete  methods. 

1.3  Related  Algorithms 

The  robotics  literature  on  path  planning  in  the  presence  of  obstacles  is  re¬ 
lated  to  the  mission  planning  problem.  Two  general  classes  of  problems  are 
addressed,  one  related  to  mobile  robots  and  the  other  related  to  industrial 
robotic  arms  (manipulators). 

Efficient  methods  for  generating  shortest  paths  for  mobile  robots  have 
been  developed.  They  are  generally  based  on  finding  the  shortest  path  in  the 
tangent  graph  (or  visibility  graph)  consisting  of  the  collison-free  tangents 
between  all  obstacles.  For  convex  polygonal  obstacles,  this  method  is  based 
on  the  premise  that  when  two  points  in  a  plane  are  not  visible  from  each  other, 
the  shortest  path  always  contains  one  or  more  vertices  of  the  obstacle.  Part  of 


a  tangent  graph  is  illustrated  in  Figure  4.  The  naobile  robot  will  move  from  M 
to  T  on  the  line  segments  composed  of  edges  of  the  polygons  and  edges  linking 
extreme  points  of  the  polygons.  This  approach  works  best  for  polygons  and 
does  not  easily  adapt  to  obstacles  having  curved  surfaces.  Although  the 
circles  may  be  approximated  by  polygonal  shapes,  a  close  approximation 
requires  many  turns.  In  the  case  of  mission  planning,  a  path  having  many 
segments  is  undesirable.  Variations  of  this  strategy  can  be  found  in  [7,  10, 
11,  12,  13]. 

Others  have  treated  this  mobile  robot  problem  as  a  control  problem. 
Shortest  paths  are  generated  by  following  the  negative  gradient  of  some  re¬ 
turn  function.  Multiple  obstacles  are  avoided  one  at  a  time  and  the  paths 
are  smooth  curves  (see  [18,  19,  26]). 


Figure  4  about  here 


Operators  of  industrial  robots  have  to  either  teach  or  program  a  manip¬ 
ulator  so  that  it  can  function  in  the  work  space  without  collisions.  In  this 
case,  the  strategies  for  path  planning  include  optimization  which  treats  the 


# 


obstacle  as  a  repulsive  force  field.  When  a  robot  arm  moves  close  to  the 
obstacle,  a  large  repulsive  force  is  generated  which  penalizes  movement  close 
to  an  obstacle.  Optimization-based  strategies  may  be  found  in  [3,  4].  Other 
strategies  include  sensor-based  motion  planning  [21,  22],  control  methods 
[15],  and  geometric  methods  [24].  An  algorithm  for  a  coordinate  measuring 
machine  which  requires  that  a  probe  be  positioned  at  various  locations  on 
a  part  may  be  found  in  [14].  These  algorithms  generally  produce  curved 
paths  and  also  involve  considerations  of  velocity,  acceleration,  joint  types, 
and  manipulator  range. 

The  problem  of  finding  shortest  paths  in  2D  with  obstacles  has  been 
addressed  by  Zheng  and  Iyengar  [25].  They  place  a  uniform  grid  over  the 
plane  and  apply  algorithms  to  the  corresponding  graph  problem.  For  the  3D 
problem,  Jun  and  Shin  [9]  divide  the  workspace  into  cells  so  that  the  problem 
becomes  discrete  with  the  obstacles  occupying  known  cells.  Using  this  idea, 
the  problem  is  to  find  a  sequence  of  neighboring  cells  from  the  origin  to 
the  destination  which  exclude  all  cells  intersecting  obstacles.  These  methods 
tend  to  find  solutions  that  only  allow  certain  angles  between  segments,  such 
as  90°  angles.  Often  such  paths  can  be  shortened,  with  fewer  segments,  if 
there  is  less  restriction  in  the  choice  of  turn  angle. 
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The  advantage  of  our  algorithm  is  that  it  accomodates  both  circular  and 
polygonal  shapes.  In  addition,  it  generates  straight  line  segments  even  for  the 
circular  shapes,  which  is  required  in  order  to  program  the  missile  waypoints. 
The  small  number  of  segments  simplifies  this  programming,  and  the  short 
paths  are  efficient  and  increase  the  probability  of  success. 

1.4  Other  Applications 

Layout  and  routing  problems  also  involve  finding  paths  which  avoid  obsta¬ 
cles.  Consider  the  problem  of  routing  fiber  optic  cable  between  major  cities. 
The  objective  is  to  determine  a  short  route  that  avoids  certain  regions  such 
as  other  cities,  airports,  military  bases,  lakes,  and  mountains.  A  similar 
problem  is  faced  by  public  utilities  (gs-s,  water,,  electric),  as  well  as  high¬ 
way  planners.  Trogneux,  Doquet,  and  Mallet  [20]  address  the  problem  of 
designing  a  high  voltage  transmission  network  which  avoids  obstacles  such 
as  airfields  and  army  bases  and  takes  environmental  obstacles  into  account. 
Microwave  communication  links  have  a  line-of-sight  requirement.  Finding 
the  placement  of  microwave  towers  between  major  cities  may  be  viewed  in 
this  way.  In  these  types  of  applications,  the  turn  angle  constraint  which  our 
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algorithm  addresses  becomes  unimportant. 

1.5  Objective  of  the  Investigation 

Advances  in  graphics  technology  has  led  to  a  new  field  of  study  known  as 
scientific  visualization  or  just  visualization.  Research  groups  in  physics,  fluid 
dynamics,  meteorology,  medicine,  molecular  biology,  and  operations  research 
are  using  visualization  techniques  (including  animation)  to  help  in  problem¬ 
solving  ([8]).  In  optimization  applications,  visualization  is  frequently  used  to 
display  solutions  to  complex  multidimensional  problems.  This  investigation 
not  only  uses  visualization  to  display  the  solutions,  but  also  uses  visualization 
to  evaluate  the  solutions.  Simply  stated,  the  objective  of  this  investigation 
is  to  develop  algorithms  which  automatically  generate  feasible  missions  for 
cruise  missiles  which  will  be  acceptable  to  mission  planners.  The  algorithm 
provides  paths  which  to  human  perception  axe  a  good  balance  betwen  short 
total  path  distance  and  a  small  number  of  segments.  Hence,  the  empirical 
analysis  results  in  a  set  of  figures  in  which  appearance  is  the  only  criteria  for 
evaluation. 
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2  CIRCULAR  THREATS 


In  this  section,  the  K  convex  sets  Ri,...,Rk  are  assumed  to  be  circles  each 
with  center  Xi,...,Xk  €  R"^  and  radii  ri,...,rA-.  The  method  described 
in  this  section  is  called  the  circumscribed  triangle  algorithm.  It  starts  at  M 
(the  missile  location)  and  constructs  line  segments  by  moving  from  a  given 
segment  point  tangentially  to  the  circumference  of  nearby  threat  regions, 
until  a  suitable  straight-line  path  from  a  segment  point  to  T  (the  target 
location)  is  found.  This  basic  strategy  is  embedded  in  a  branch-and-bound 
framework,  so  that  a  short  path  of  this  restricted  type  is  obtained. 

In  this  algorithm  we  allow  the  missile  to  intersect  a  point  on  the  circum¬ 
ference  of  a  threat  region  but  do  not  allow  the  missile  to  fly  through  any 
interior  point  of  a  threat  region.  Given  segments  and  [y;,yi+i],  the 

turn  angle  at  segment  point  Yi  is  defined  to  be  the  complement  of  angle 

Let  d{E,  F)  denote  the  Euclidean  distance  between  points  E  and  F.  The 
basic  step  in  proceeding  from  segment  point  Yi  tangentially  to  a  (nearby) 
threat  region  Rk  =  {s  :  d{Xk,s)  <  r*}  involves  constructing  the  circum¬ 
scribing  triangle  on  Rk  having  one  vertex  at  Yi  (see  Figure  5). 


Figure  5  about  here 


# 


The  tangential  points  in  Figure  5  are  Ti  ,  T2,  and  T3  and  the  vertex  points 
are  Yi,  Fi,  and  V2.  Sides  YiVi  and  YV2  are  the  equal  sides  of  the  (isosceles) 
circumscribing  triangle. 

The  tangential  points  in  Figure  6  are  Ti  ,  T2,  and  T3  and  the  vertex  points 

are  Yi,  Vi,  and  V2.  Sides  Y{Vi  and  YV2  are  the  equal  sides  of  the  (isosceles) 

circumscribing  triangle.  Tangential  point  Tz  is  easily  computed  since  it  lies 

on  the  line  through  Y,  and  Xk  with  d(Xk,  T3)  =  r^.  Auxiliary  point  Z  formed 

by  dropping  a  normal  to  YjTz  from  Ty  (or  r2)  can  be  used  to  obtain  Ti  and 

Vi  (or  T2  and  V2).  Note  that  AYiZTy,  AYiTzVy,  AYTyXk,  and  AZTyXu  are 

all  similar  triangles. 

Since  AZTyXk  -  AYiTyXu, 

XkZ  XkTy  d{Xk,Z)  _  Tk  _ 

XtT,  “  XtYi  ”  r*  d{Yi,X,) 

so  that 

=  (rtfld(Yi,Xt)  . 
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Also  since  AZTiXk  ~  AYiTiXk, 


ZT,  TiYi  d{TuZ)  _  \ld(Yi,XtY-(nf 
Xtf;  ~  XtYi  “  n  “  il.Yi,Xt) 


SO  that 

d{T^,Z)  =  rky/d{Yi,Xky  -  {rk)yd{Yi,Xk)  . 


Since  AYT3V1  ~  AYiZTu 


ViTz  TiZ  d{Ts,Vi)  _  djTi^Z) 

YiTz  ~  YiZ  “  d{YuXH)  +  rk  d{Y,Xk)  -  d{Xk,Z) 


SO  that 


d{T^,V^) 


{d{Yi,Xk)  +  n)n-^d{Yi,Xky  -  [rkf 
d{Yi,Xk){d{Xi,  Xk)-d{Xk,Z)) 
{d{Yi,Xk)  +  n)n^Jd{Yi,Xky-{rky 
d{Yi,Xk){d{Yi,Xk)  -  {rkyid{Yi,Xk)) 
d{Yi,Xk)n)  +  {rky 
yld{Yi,Xky-{rky  ■ 


Figure  6  about  here 


Thus  points  Vi  and  V2  can  be  readily  obtained  by  moving  a  distance 
d{T3,Vi)  along  the  two  normals  to  Y1T3  at  point  T3.  The  points  Ti  and 
can  be  obtained  by  moving  a  distance  yJd{Y{,Xky  —  (rjt)^  from  Yi  toward  Vi 
and  Vz,  respectively,  or  by  computing  the  nearest  points  to  Xk  lying  on  the 
lines  through  Y,  and  Vi  or  V2. 

We  intend  to  make  use  of  these  tangential  lines  to  construct  the  next 
segment,  if  possible.  Before  using  points  on  either  tangential  line,  we  first 
make  sure  that  the  respective  turn  angles  at  Yi  would  be  within  the  turn 
angle  limit.” 

The  points  Vi  and  V2  have  a  special  property.  Each  is  the  first  point  on  the 
tangential  line  from  Yi  through  itself  which  has  a  complete  180°  field  of  view 
(ignoring  turn  angle  constraints)  of  the  side  of  Rk  opposite  Yi.  By  contrast, 
tangent  points  Tj  and  T2  each  have  less  than  a  90°  field  of  view  (ignoring 
turn  angle  constraints)  of  that  side  (see  Figure  7).  With  this  in  mind  Vi  (or 
V2)  would  seem  to  be  a  good  choice  for  the  end  point  of  a  next  segment  in 
our  path  construction.  We  actually  attempt  to  use  a  point  between  Ti  and 
Vi  (or  between  T2  and  V2)  for  the  next  end  point.  In  general,  such  a  point 
can  be  represented  as  a  convex  combination  Yi+i  =  (1  —  a)Ti  +  aVi  (or 
Yi+i  =  (1  —  Q:)r2  +  C(V2  ),  for  some  0  <  a  <  1.  In  order  for  such  a  point 


Yi+i  to  be  usable,  we  must  insure  that  the  segment  is  clear.  (A  line 

segment  is  defined  to  be  clear  if  it  does  not  intersect  the  interior  of  any  of 
the  threat  regions.) 


Figure  7  about  here 


We  first  search  for  a  minimum  value  of  a  with  0  <  or  <  1  such  that 
[Yi,  (1  -  Q:)ri  +  aVi]  (or  [Fi,  (1  -  «)22  +  ocVt]  )  is  clear  and  [(1  -  a)Ti  +  aVi,T] 
(or  [(1  —  Q;)T2  +  aV2,T]  )  is  clear  without  exceeding  the  turn  angle  limit  at 
(1  -  a)Ti  +  aVi  (or  (1  -  a)r2  +  aV2  ).  If  such  a  point  is  found  we  have  a 
feasible  completion  by  adding  the  two  segments  [Yi,Yi+i  =  (1  —  a)Ti  +  aVi] 
and  [F;-+x  =  (1  -  a)Ti  +  aVi, ^^+2  =  T]  (or  [F;-,  F;-+i  =  (1  -  a)T2  +  ay2]  and 
[F:-+i=(l-a)r2  +  r.V2,F^-+2  =  T]). 

If  such  a  point  leading  to  a  feasible  completion  cannot  be  found,  it  would 
seem  natural  to  use  Vi  (or  V2)  by  default,  due  to  its  wide  field  of  view. 
However,  in  practice  this  leads  to  longer  paths  and  we  find  it  better  to  use  a 
point  between  2i  and  Vi  (or  between  T2  and  V2)  for  the  next  end  point.  We 
have  parameterized  this  by  defining  a  backoff  percentage  ^  which  specifies 


that,  in  this  case,  Vi+i  is  to  be  chosen  as  the  convex  combination  (1  -  ^)Tx  + 
(or  (1  -  P)T2  +  PV2  ).  Experimentally,  we  have  found  /?  =  |  to  be  a  good 
choice. 

The  branch-and-bound  algorithm  we  have  developed  based  on  this  cir¬ 
cumscribed  triangle  approach  is  given  as  follows: 


Algorithm  :  circumscribe 
Inputs  : 


M 

T 

K 

Xi,...,Ak 

ri,...,rK 

Z 


-  missile  location 

-  target  location 

-  number  of  threats 

-  center  points  of  the  threats 

-  radii  of  the  threats 

-  turn  angle  limit  (degrees) 

-  fall  back  parameter  (0  <  /?  <  1) 


Outputs  : 

D 

S 


-  length  of  best  constructed  solution  path 

-  number  of  segments  in  best  constructed  solution  path 


Li, Ls 

-  segments  of  best  constructed  solution  path 

Constructs  : 

B 

-  pointer  to  last  point  in  best  constructed  solution  path 

Q 

-  branch-and-bound  queue 

diE,F) 

-  distance  from  E  to  F 

L{E,F) 

-  line  segment  from  E  to  F 

R{C,r) 

-  threat  region  with  center  C  and  radius  r 

S{C,r) 

-  interior  of  threat  region  with  center  C  and  radius  r 

c 

I 

-  number  of  points  constructed 

Yi 

-  constructed  point 

Bi 

-  pointer  to  point  previous  to  Yi 

Di 

-  path  distance  prior  to  Yi 

Si 

-  path  segments  prior  to  Yi 

Ei 

-  set  of  threat  points  eligible  for  moving  away  from  point  5^ 

-  constructed  tangent  points  of  circumscribled  triangle 

Vr,V2 

-  constructed  vertices  of  circumscribled  triangle 

begin 

D  *—  oo;  S  * —  0;  5  <—  0; 
if  L{M,  T)  n  S{Xi,  r.)  =  $  for  i  =  I..K  then 
D^d(M,T);  S^l]  Li<-[M,r]; 


begin 

/^l;  Si^O;  Di^O;  5i  ^  0;  Ei  *- {Ru-,Rk} 

put  {Yi,Bi,Di,Si,Ei)  on  Q-, 
while  (5  ^  $  do 
begin 

remove  some  (yj,  Bi,  D{,  Si,  Ei)  from  Q; 
b  *—  Bi] 

select  some  Rj  from  Ei;  Ei  ^  Ei\  Rj; 
if  Ei  ^  $  then  put  {Yi,Bi,Di,  Si,Ei)  back  on  Q; 
form  the  triangle  circumscribing  Rj  having 

its  tangent  points  Ti  €  YVi,T2  €  YiV2,Tz  €  ViV^; 
for  k  =  1  to  2  do 
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begin 

iim-YbYiTk<Z  and 

L{Yi,Tk)^S{Xi,ri)  =  ^  for  /  =  1../^  and 
A  +  d{Yi,Tk)  +  d{Tk,T)  <  D  then 
begin 

search  for  a  point  U  on  L{Tk,  T4)  nearest  to  Tk  satisfying: 
180  -  YiUT  <  Z  and 

L{Yi,U)r]S{Xi,ri)  =  ^  for  I  =  1..K  and 
L{U,T)r\SiXi,ri)  =  ^  for  l=l..K-, 
if  such  a  point  U  exists  then 
begin 

if  A  +  d{Yi,  U)  +  d{U,  T)<D  then 
begin 

Yi^U\  Bi^i]  Di  ^  Di  +  d{Yi,Uy,  Si^Si  +  Y, 

m  *—  I\  7  <—  /  +  1; 

y/^T;  Di  +  d{U,Ty,  Si  ^  Si +  2; 


A-22 


D  ^Di-  S  *- Sr,  B^r, 
end; 
else 
begin 

if  Ei  ^  $  then 
begin 

U  ^  (I- l3)Tk  +  0Vk-, 

if  L{Yi,U)nS{Xi,ri)  =  ^  for  /  =  1../^  then 
begin 

/<—/  +  !; 

Yi  V,  Dj  4-  Di-Vd{Yi,  U);  Si  <-  5,-+l;  Ei 

put  (Yi,  Bi,  Di,  Si,  Ei)  on  Q-, 
end; 
end; 


end; 

end; 

end; 


end; 

end; 
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while  k  ^  0  do  j  <—  Bk\  Li  <—  [Yj,  Yk]',  i  •<—  i  —  1;  <—  i; 

end; 

Figure  8  shows  an  optimal  path  through  a  region  containing  three  threats 
generated  by  the  circumscribed  triangle  algorithm.  Figure  9  shows  additional 
partial  paths  produced  by  the  branch-and-bound  process  contained  in  the  al¬ 
gorithm.  The  actual  implementation,  whose  results  are  reported  in  Appendix 
A,  actually  makes  two  applications  of  the  algorithm,  the  second  of  which  has 
the  roles  of  missile  and  target  reversed.  Of  course,  any  path  obtained  with 
missile  and  target  reversed  would  itself  need  to  be  reversed.  The  shortest  of 
the  two  paths  generated  is  selected  as  the  best  path. 


Figures  8  and  9  about  here 


3  CIRCULAR  THREATS  AND  A 
PROBABILITY  CONSTRAINT 

In  this  section,  a  threat  is  represented  by  a  circle  with  a  given  center  X  and 
a  given  radius  r.  Let  Q  denote  the  point  along  the  mission  plan  at  which  the 
missile  is  nearest  X.  If  \\Q  —  ^||  >  t*,  then  the  threat  at  A  is  too  far  away  to 
harm  the  missile  and  the  probability  of  a  successful  traversal  is  defined  to  be 
1.  Otherwise,  some  function  p{Q,  X,r)  is  used  to  determine  the  probability 
of  a  successful  traversal.  When  there  are  multiple  threats,  the  probability  of 
a  successful  traversal  is  the  product  of  the  probabilities  along  each  segment, 
teiking  into  account  all  of  the  threats. 

The  problem  of  finding  a  feasible  mission  plan  which  incorporates  a  prob¬ 
ability  constraint,  can  be  defined  mathematically  as  follows: 

Given  a  target  probability  P,  two  points  M  and  T  ^  ,  and  K 

circles,  with  centers  Ai, ...,  Aa'  €  and  radii  ri,  ...,rA', 
find  a  set  of  line  segments 

with  Yi  =  M  and  Ys+i  =  T,  such  that  the  probability  of  a 
successful  traversal  is  at  least  P. 

The  quality  of  a  mission  plan  is  measured  by  the  length  of  the  path  and  the 


number  of  path  segments.  The  best  path  is  the  straight  line  from  M  to  T.  If 
this  is  not  feasible  (i.e.  the  probability  of  a  succcessful  traversed  is  less  than 
P),  then  we  seek  a  short  path  having  only  a  few  segments. 

The  probabilistic  incursion  triangle  algorithm  allows  the  missile  to  fly 
through  the  interior  of  a  threat  region  and  employs  a  probability  model  to 
compute  the  probability  of  survival  for  a  path  from  the  missile  location  M 
to  the  target  location  T  which  consists  of  line  segments  which  may  either  be 
tangent  to  a  circular  threat  region  or  cut  through  the  region  (an  incursion). 

The  algorithm  operates  within  the  same  branch-and-bound  algorithmic 
structure  used  in  the  circumscribed  triangle  algorithm.  The  major  differences 
are  (1)  when  attempting  to  construct  a  new  segment  by  approaching  a  threat 
region  from  a  previous  segment  endpoint,  several  equilateral  triangles  passing 
through  the  threat  region  in  addition  to  the  circumscribed  triangle  may  be 
considered,  and  for  each  such  triangle,  potentially  two  points  are  added  to 
the  branch-and-bound  queue,  and  (2)  in  searching  for  a  final  segment  linking 
a  previous  segment  endpoint  to  the  target,  that  segment  may  intersect  threat 
regions  (with  a  consequent  reduction  in  survival  probability). 

For  threat  region  k  with  center  Xk  and  radius  r^,  the  probability  model 
assumes  that  the  probability  of  survival  of  a  line  segment  approaching  a 


threat  region  is  a  function  only  of  the  minimum  distance  from  the  segment 
to  Xk-  Further,  the  model  assumes  independence  so  that  the  overall  survival 
probability  for  a  path  is  the  product  of  the  survival  probabilities  for  each 
segment  in  the  path. 

For  developmental  purposes  we  have  hypothesized  a  simple  probability 
function  which  is  a  quadratic  of  the  form 

Pr{d)  =  aS  +  c  , 

where  for  threat  region 

Pr(rjt)  =  .99  and  Pr(0)  =  .50  . 

Figure  10  illustrates  a  circumscribed  triangle  YiV\V2  an  incursion 
triangle  YiQiQ2  for  a  common  threat  region  Rk-  In  the  circumscribed  triangle 
algorithm,  with  respect  to  the  segments  K'Vi  and  3^14,  the  line  segments 
T\V\  and  T2V2  are  searched  first  for  a  clear  completion  segment  to  the  target, 
before  the  points  Vi  and  V2  are  added  to  the  branch-and-bound  queue. 

This  also  occurs  in  the  probabilistic  incursion  triangle  algorithm,  but  the 
completions  can  intersect  threat  regions  if  the  total  survival  probability  is  at 
least  P  and  the  segments  YiVi  and  YiV2  will  each  have  survival  probability 
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.99.  Also,  in  the  probabilistic  incursion  triangle  algorithm,  with  respect  to 
the  segments  YiQi  and  Y{Q2,  the  line  segments  PiQi  and  P2Q2  are  searched 
first  for  a  completion  to  the  target  with  total  survival  probability  at  least  P, 
before  the  points  Vi  and  Vi  are  added  to  the  branch-and-bound  queue. 


Figure  10  about  here 


The  number  of  incursion  triangles  considered  in  approaching  a  threat 
region  from  a  previous  segment  endpoint  depends  on  the  accumulated  number 
of  segments  and  an  adjustable  parameter  S  which  is  roughly  interpreted  as 
the  expected  number  of  segments  in  an  optimal  path.  In  the  testing  described 
in  Appendix  B,  S  was  set  at  3.  If  Si  is  the  number  of  segments  in  a  partial 
path  from  the  missile  location  to  segment  endpoint  Yi  and  Pi  is  the  survival 
probability  of  this  partial  path,  incursion  triangles  with  individual  segment 
survival  probabilities  of 

M{p/Pi) ,  .m . ) 


will  be  considered,  unless  these  probabilities  are  very  close  to  .99. 


The  twenty  test  problems  presented  in  Appendix  A  were  solved  using 
P  =  90%,  P  =  75%,  and  P  =  60%,  and  the  results  are  summarized  in  Table 
1.  The  deterministic  algorithm  is  the  circumscribed  triangle  algorithm.  Note 
that  setting  P  =  60%  subsantially  reduces  the  number  of  waypoints,  but  only 
reduces  the  total  distance  by  4%.  The  paths  obtained  by  the  deterministic 
algorithm  and  the  new  algorithm  with  P  =  60%  may  be  found  in  Appendix  B. 
For  most  problems,  the  two  paths  are  similar,  with  the  new  algorithm  taking 
short  cuts  through  one  or  more  threats.  Problem  7  illustrates  a  case  where 
the  new  algorithm  selects  a  completely  different  path.  For  this  problem, 
the  deterministic  algorithm  produces  a  path  in  which  all  threats  are  above 
the  path,  while  the  new  algorithm  uses  a  mixed  strategy.  The  deterministic 
strategy  is  unable  to  use  a  single  straight  line  for  any  of  these  test  problems, 
whereas  the  probabilistic  algorithm  finds  this  to  be  best  for  several  of  them. 


Table  1  about  here 


4  POLYGONAL  THREATS 


In  this  section,  we  extend  the  ideas  presented  in  Section  2  to  handle  the  case 
of  threat  regions  modelled  as  convex  hulls.  The  algorithm  described  below  is 
fundamentally  the  same  as  the  algorithm  described  in  Section  2.  The  modi¬ 
fied  algorithm  is  given  in  Appendix  C.  The  original  algorithm  required  that 
all  threat  regions  be  modelled  as  circular  regions.  The  major  modification 
we  have  made  also  allows  a  threat  region  to  be  modelled  as  the  convex  hull 
of  a  given  set  of  points. 

The  original  and  modified  algorithms  start  at  M  (the  missile  location)  and 
construct  line  segments  by  moving  from  a  given  segment  point  tangentially 
to  the  boundary  of  a  nearby  threat  region.  This  process  continues  until  a: 
suitable  straight-line  path  from  a  segment  point  to  T  (the  target  location)  is 
found.  This  basic  strategy  is  embedded  in  a  branch-and-bound  framework, 
so  that  a  short  path  of  this  restricted  type  is  obtained. 

The  missile  is  allowed  to  intersect  a  point  on  the  boundary  of  a  threat 
region  but  is  not  allowed  to  fly  through  any  interior  point  of  a  threat  region. 
Further,  a  turn  angle  limit  is  enforced  at  each  segment  point,  since  a  missile 
is  unlikely  to  be  making  a  sharp  turn  at  such  a  point.  Given  segments 


[Fi-ijli]  and  [K-,5^+i],  the  turn  angle  at  segment  point  Yi  is  defined  to  be 
180°  —  lYi-iYiYi^i. 

Let  d{E,  F)  denote  the  Euclidean  distance  between  points  E  and  F.  The 
basic  step  in  proceeding  from  segment  point  Yi  tangentially  to  a  (nearby) 
threat  region  Rk  involves  constructing  the  circumscribing  triangle  for  Rk 
which  has  one  vertex  at  Yi.  Previously,  a  procedure  was  described  for  ob¬ 
taining  the  circumscribing  triangle  used  when  moving  from  a  segment  point 
Yi  tangentially  past  a  circular  threat  region  defined  by  a  center  point  Xk  and 
radius  r^. 

The  modified  algorithm  can  also  obtain  the  circumscribing  triangle  to  be 
used  in  moving  from  a  segment  point  Yi  past  a  threat  region  Rk  defined  as 
the  convex  hull  of  a  set  of  generator  points  {C?i, . . . ,  Gj}  (see  Figure  11). 


Figure  11  about  here 


The  generator  points  defining  the  circumscribed  triangle  in  Figure  11  are 
the  side  points  Ti  ,  T^,  and  Ts  (which  play  a  role  similar  to  that  of  the 


tangential  points  in  the  case  of  a  circular  region)  and  the  vertex  points  are 
Yi,  Vi,  and  V^.  Sides  YiVi  and  YV2  are  the  equal  sides  of  the  (isosceles) 
circumscribing  triangle.  Points  Ti  and  are  the  generator  points  associated 
with  the  largest  angle  on  distinct  generator  points  having  vertex  Yi.  This 
angle  is  actually  determined  by  computing 

{Ga  -  Y)  •  (Gi,  -  Yi) 

I  (G,  -  Yi)  II  (Gi  -  «)  I  ■ 

The  generator  point  Tz  is  associated  with  the  largest  inner  product  formed 
from  generator  points  and  (normal  vector) 

(Ti-Yi)  {Tz-Yi)  \ 

l(ri-y:ori(r2-T;-)i;  ' 

Note  that  the  line  from  Yi  in  the  direction  given  by  the  above  normal  vector 
bisects  the  angle  TiYiTz  and  determines  the  auxiliary  point  Z,  which  in  turn 
bisects  the  segment  V1V2.  The  maximum  inner  product  is  computed  as 

77  =  max  Ga  • 
a 

Points  Vi  and  V2  are  then  given  by 


(Tr-Y)  .  {T2-Y)  \ 

\{T^-Y)\^  \{T2-Y)\)  ■ 


Vi 


t;  + 


ai 


(  in-Y)  \ 
\\{n-Y)\) 
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where  ai  and  are  determined  from  the  solution  of  the  equations: 


V 


Yi  +  ai 

Yi  +  02 


/'jJWilL'l'l 


(T.-y<)  ,  (T2-Yi)  \ 

l(r.-y)l  I(r2-v;)i; 
(r.-y.)  ,  (r^-y)  'i 
l(T.-y)r  i(rj-K)i;  ’ 


which  place  Vi  and  V2  on  the  line  normal  to  YiZ  at  Z.  Note  that  Tz  is  also 
on  that  normal  line. 

Let  W\  (or  W2)  be  a  fixed  point  along  the  line  segment  starting  at  T\  (or 
T2)  through  V-i  (or  V2).  The  algorithm  searches  for  a  minimum  value  of  a 
with  0  <  a  <  1  such  that  [3^-,  (1  — Qr)ri+Q;lLi]  (or  [3^-,  (1  — a)T2  +  o:H'2]  )  does 
not  intersect  the  interior  of  any  of  the  threat  regions,  [(1  —  Q;)ri  +  aWi^T\ 
(or  [(1  —  Oi)T2  +  0^2)  T]  )  does  not  intersect  the  interior  of  any  of  the  threat 
regions,  and  the  turn  angle  limit  was  not  exceeded  at  (1  —  a)ri  +  aW\  (or 
(1  —  a)r2  +  aW2  ).  In  the  original  algorithm,  W\  was  Vi  and  W2  was  V2" 
The  modified  algorithm  performs  this  seaxch  for  both  a  circular  threat  region 


and  a  convex  hull  threat  region,  and  if  the  line  segment  from  Yi  through  Ti 
(or  T2)  to  the  boundary  does  not  interesect  any  of  the  threat  regions,  then 
Wi  (or  W2)  is  taken  to  be  that  boundary  point  instead.  If  such  a  suitable 
point  is  found,  a  feasible  completion  is  formed  by  adding  the  two  segments 
[Yi,Yi+i  =  (1  -  Ci)Ti  +  aWi]  and  [K'+i  =  (1  -  <y)Ti  +  aWi,yf+2  =  (or 
[li,  yj+i  =  (1  —  <y)T2  +  0LW2]  and  [K+i  =  (1  ~  0‘)T2  +  0:^2, 5^+2  =  2^])  to  the 
line  segments  already  constructed  leading  from  M  to  K'. 

If  such  a  point  leading  to  a  feasible  completion  cannot  be  found,  a  point 
between  T\  and  Vi  (or  between  T2  and  V2)  is  used  for  the  next  end  point. 
In  the  original  algorithm,  a  backoff  percentage  ^0  =  50%  was  used  so  that 
the  next  segment  point  Yi+i  is  the  convex  combination  {1  -  ^)Ti  +  fiVi 
(or  (1  -  0)T2  +  I3V2  ).  In  the  modified  algorithm,  ^  is  computed  so  that 
(1  —  P)Ti  +  fSVx  (or  (1  -  P)T2  +  ^¥2  )  is  the  first  point  from  Ti  toward  Vi 
(or  from  T2  toward  V2)  such  that  a  line  at  that  point  with  the  limiting  turn 
angle  intersects  the  boundary  of  the  threat  region,  but  not  the  interior.  (For 
a  convex  hull  threat  region,  a  degenerate  case  is  possible  in  which  T\  =V\ 
or  T2  =  V2  and  a  distinct  intermediate  point  is  neither  possible  nor  needed.) 

As  with  the  original  implementation,  the  algorithm  is  applied  twice,  the 
second  time  with  the  missile  and  target  locations  reversed,  and  the  better  of 


the  two  paths  generated  is  taken  as  the  recommended  path. 

The  real  test  for  this  algorithm  is  whether  or  not  it  develops  paths  which 
look  reasonable  when  displayed  on  a  two-dimensional  drawing.  The  paths 
for  twenty  problems  are  displayed  in  Appendix  D.  These  paths  not  only 
appear  reasonable,  but  also  appear  to  be  of  minimum  length.  In  addition, 
the  solution  time  on  a  Dec  5000  workstation  was  less  than  one  minute  for 
each  of  these  test  problems. 

5  SUMMARY  AND  CONCLUSIONS 

This  manuscript  presents  a  new  algorithm  to  automatically  generate  a  path 
composed  of  line  segments  from  a  given  origin  (missile  location)  to  a  given 
destination  (target  location)  which  avoids  obstacles  (threat  regions)  modelled 
as  circles.  The  algorithm  exploits  the  geometry  of  the  problem  description 
and  mimics  the  trial-and-error  strategy  a  human  mission  planner  might  wish 
to  apply  to  this  graphical  problem.  The  basic  idea  is  that  short  paths  which 
avoid  circular  obstacles  should  coincide  with  the  edges  of  a  circumscribing 
triangle  defined  by  a  current  point  and  a  circle  between  the  current  point  and 
the  target.  A  heuristic  strategy  is  applied  to  determine  the  stopping  point 


along  each  of  these  two  sides.  Ideally  the  stopping  point  will  have  a  clear 
path  to  the  target,  in  which  case  a  path  has  been  discovered.  If  this  is  not 
the  case,  then  this  is  repeated  with  each  of  the  previous  stopping  points.  At 
each  step  there  are  two  edges  which  can  be  followed  to  obtain  a  path  around 
an  obstacle.  All  such  paths  are  examined  and  the  shortest  is  retained.  The 
algorithm  reverses  the  roles  of  the  missile  and  the  target  and  the  shortest 
path  from  each  application  is  the  path  reported  by  the  algorithm. 

Two  extensions  of  this  basic  algorithm  are  also  presented  in  this  manuscript. 
One  extension  assigns  a  probability  of  a  safe  traversal  for  each  leg  of  a  path 
from  the  origin  to  the  destination.  Legs  are  allowed  to  pass  through  the 
threat  circles  with  a  corresponding  reduction  in  the  probability  of  a  safe 
traversal  along  this  leg.  The  algorithm  is  extended  to  produce  short  paths 
which  satisfy  a  given  probability  restriction.  The  second  extension  is  for 
threat  regions  which  are  convex  polygons.  The  mission  planning  problem  as 
described  in  this  manuscript  has  no  objective  function  and  can  only  be  judged 
qualitatively  by  examination  of  a  graphical  display  of  the  problem  and  solu¬ 
tion  path.  We  believe  that  the  sixty  displays  presented  in  this  manuscript 
provide  a  convincing  case  that  our  algorithms  provide  good  solutions. 
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a.  Shortest  Path 


Figure  2: 


A  Mission  with  a  Single  Threat 
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Figure  3:  A  Mission  with  Two  Threats 
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Figure  9:  Partial  Paths  Considered  by  the  Circumscribed  Triangle 
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Table  1.  Comparison  of  Probabilities  for  the  Mission  Planning  Problem 


Prob 

Deterministic 

Probabilistic  Algorithm 

Algorithm 

90% 

75% 

60% 

Dist 

Segmts 

Dist 

Segmts 

Dist 

Segmts 

Dist 

Segmts 

1 

957 

3 

955 

3 

952 

3 

949 

3 

2 

853 

3 

849 

3 

842 

3 

834 

3 

3 

853 

4 

851 

3 
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3 
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3 

4 

945 

3 

939 

3 
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3 

905 

3 

5 

795 

2 
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2 

789 

2 

783 

2 

6 

878 

4 
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4 

868 

4 
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4 

7 

1042 

3 

1040 

3 

1031 

3 

1012 

3 

8 

805 

2 
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2 
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1 
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1 

9 

782 

2 

782 

2 

780 

2 

780 

1 

10 

1049 

3 

1044 

3 

1031 

3 

1012 

3 

11 

847 

3 

845 

3 

838 
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835 

2 

12 

820 
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819 
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817 
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815 

1 

13 

°810 
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809 

2 

807 

2 
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1 

14 

945 

3 

942 

3 

861 

3 

809 

2 

15 

966 

3 

962 

3 

950 

3 

780 

2 

16 

837 

3 

837 

3 

835 

3 

835 

1 

17 

865 

3 

857 

3 

848 

3 

825 

3 

18 

1009 

3 

1005 

3 
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3 

977 

3 

19 
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4 

895 

3 
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3 
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3 

20 

862 

2 

851 

3 

837 

2 

820 

3 

Totals 

17820 

57 

17757 

56 

17534 

54 

17099 

47 

Scaled 

1.0 

0.982 

0.984 

0.947 

0.960 

0.825 

APPENDIX  A 


SOLUTIONS  OBTAINED  USING  THE 
CIRCUMSCRIBED  TRIANGLE  ALGORITHM 

Using  a  grid  size  of  780  x  700,  we  randomly  generated  five  groups  of  test 
problems.  Each  group  had  five  circular  threats  with  radii  randomly  generated 
on  the  interval  [50, 150].  The  four  problems  in  each  group  had  identical  threat 
sites,  but  different  missile  and  target  locations  with  the  missiles  located  on 
the  left  boundary  and  targets  located  on  the  right  boundary.  The  missile 
and  target  locations  were  randomly  selected  and  were  different  for  each  of 
the  twenty  test  problems.  The  input  for  each  problem  consists  of  seven  points 
(the  missile  location,  the  target  location,  five  threat  centers)  and  five  radii. 

A  graphical  system  was  developed  using  Tcl  and  Tk  (see  Welch  [23]). 
Tcl  stands  for  Tool  Command  Language  and  Tk  is  a  toolkit  for  window  pro¬ 
gramming.  The  system  allowed  us  to  display  the  strike  mission  graphically. 
The  paths  obtained  by  the  circumscribed  triangle  algorithm  may  be  found 
in  Figures  A.l  through  A.5. 


Figure  A.l:  Display  of  Problems  1,  2,  3,  and  4 


Problem  5 


Problem  6 


Figure  A. 2:  Display  of  Problems  5,  6,  7,  and  8 
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Figure  A.3:  Display  of  Problems  9,  10,  11,  and  12 
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APPENDIX  B 

SOLUTIONS  OBTAINED  USING  THE 
PROBABILISTIC  ALGORITHM 

This  appendix  presents  a  graphical  display  of  the  twenty  test  problems  along 
with  the  solutions  obtained  by  the  deterministic  algorithm  and  the  proba¬ 
bilistic  algorithm.  The  black  line  is  the  solution  obtained  with  P  set  to  60% 
and  the  gray  line  is  the  mission  obtained  by  the  deterministic  algorithm. 
The  plots  were  obtained  using  Tcl  and  Tk  Solver  [23],  which  is  part  of  our 
experimental  computational  package. 


Problem  1 


Problem  2 


Problem  5 


Problem  6 


Figure  B.2:  Display  of  Problems  5,  6,  7,  and  8 
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Problem  13 


Problem  14 
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Figure  B.5;  Display  of  Problems  17,  18,  19,  and  20 


APPENDIX  C 

Algorithm  for  Threats  Modeled  as  Circles  and  Polygons 

The  branch-and-bound  algorithm  we  have  developed  based  on  the  circum¬ 
scribed  triangle  approach  is  as  follows: 


Algorithm: 

CIRCUMSCRIBE 


Inputs: 


M 

T 

K 

Ri,...,Rk 

Z 


-  missle  location 

-  target  location 

-  number  of  threats 

-  threat  regions 

-  boundary  of  the  rectangular  planning  area 

-  turn  angle  limit  (degrees) 


Outputs: 

D  -  length  of  best  constructed  solution  path 

S  -  number  of  segments  in  best  constructed  solution  path 

L\,...,Ls  -  segments  of  best  constructed  solution  path 


Constructs: 

B 

Q 

d{E,F) 

L{E,F) 

H{E,F) 

i?;. 

I 

Yi 

Bi 

Di 

Si 

Ei 

Ti,T2,T3 

Vl,V2 

V,W 

Pk 


-  pointer  to  last  point  in  best  constructed  solution  path 


-  branch-and-bound  queue 

-  distance  from  E  to  F 

-  line  segment  from  E  to  F 

-  halfline  from  E  through  F 

-  interior  of  threat  region  Ri 

-  number  of  points  constructed 

-  constructed  point 

-  pointer  to  point  previous  to  Y) 

-  path  distance  prior  to  Yj 

-  path  segments  prior  to  Y) 

-  set  of  threat  regions  eligible  for  moving  away  from  point 

-  constructed  tangent  points  of  circumscribled  triangle 

-  constructed  vertices  of  circumscribled  triangle 


-  vertices 

-  fall  back  fraction  {Q  <  Pk  <  1)  for  segment  TkVk 


Yi 
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begin 

D  •*—  oo;  5^0;  B  «—  0; 
if  L(Af,T)n  =  $  for  i  =  then 
D  -  d(M,Ty,  5-1;  Li  -  [Af.T]; 
else 

/  4-  1;  Yi  ^  M\  Bi  ^  0;  ^0;  5i  —  0;  ^  {Bi , Bk}; 

put  (Vi fBijDijSi,  El)  on  Q\ 
while  Q  5<i  ^  do 
begin 

remove  some  (y»,  Bi,  Di,  Si^Ei)  from  Q; 
b^Bi; 

select  some  Rj  from  Ei\  Ei  *—  Ei\  Bj-; 

if  Ei  5^  $  then  put  {Yi,  Bi,  Di,  Si,  Ei)  back  on  Q; 

form  the  triangle  YiViV^  circumscribing  Rj  having 

its  tangent  points  Tj  6  YiVi,T2  €  YiV2,T^  €  ^1^2) 

for  Jb  =:  1  to  2  do 

begin 

if  180  -  YbYiTk  <  Z  and 

L{Yi,Tk)  n  Bj)  =  ^  for  1  =  l-.B  and 
Bi  +  d(yi,Tfc)  +  d(Tfc,T)<  D  then 
begin 

W  ^  H{E,F)r\n 

ifL(yi,Tfc)nBj)  =  $  for  /  =  1..B  then 

V 
else 

V  -  n 

endif 

search  for  a  point  U  on  L(Tfc,  Y)  nearest  to  satisfying: 

180  -  YiUT  <  ^  and 
LiYi,  LT)  n  Bj)  =  $  for  /  =  1..B  and 
Llu,T)nR[)-^  for  l  =  l..K; 
if  such  a  point  U  exists  then 
begin 

ifBi  +  d(yi,Cf)  +  d(C/,T)  <  D  then 
begin 

/  ^  /+  1; 

Yj  ^U',  Bi  ^  V,  B/  -  £>i  +  d(yi,  t^);  Si  ^  Si  ^  1; 
m  -  7  +  1; 

Yi^-T]  Bf-m;  Di  ^  Dm  +  diU,Ty,  5/  ^  5i  +  2; 

D  ^Df,  s  ^Sr,  B  -  /; 

end 

else 

begin 

if  Ei  ^  ^  then 
begin 

compute  such  that  C7  ♦—  (1  -  /?)Tfe  +  ^Yk  and  180  —  yiCfTi 

ifL(Yi,C/)nB{)  =  ^  for  1  =  1.. if  then 

begin 

YrJ(/;\  -  i;  Di^Di-h  d(Yi.D);  Si  ^  Si  +  1;  Ei 
put  (y/,Br,B/,5/,E/)  on  Q; 
end 
end 
end 
end 
end 
end 
end 

i  -  S;  ib  ^  B; 

while  1:  5^  0  do  j  ♦—  Bk;  Li  ^  [y^t  ^fc]*»  *  ♦—  »  —  1;  k  j\ 
end 


=  Z; 

*-Ei; 
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APPENDIX  D 


SOLUTIONS  OBTAINED  USING  THE 

EXTENDED  CIRCUMSCRIBED 
TRIANGLE  ALGORITHM 

Using  a  grid  size  of  780  x  700,  five  groups  of  test  problems  were  randomly 
generated.  Each  group  has  five  threat  regions  composed  of  a  combination  of 
convex  sets  and  circles.  The  four  problems  in  each  group  have  the  same  threat 

c 

regions,  but  different  missile  and  target  locations.  The  missile  is  located  on 
the  left  boundary  and  the  target  is  located  on  the  right  boundary.  The  origin 
and  destination  for  each  path  were  selected  at  random  and  are  different  for 
each  of  the  twenty  test  problems. 

This  appendix  presents  a  graphical  display  of  the  twenty  test  problems 
along  with  the  solutions  obtained  by  the  algorithm.  The  plots  were  obtained 
using  Tcl  and  Tk  Solver  [23],  which  is  part  of  our  experimental  computa¬ 
tional  package. 
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Problem  5 


Problem  6 


Figure  D.2:  Display  of  Problems  5,  6,  7,  and  8 
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Problem  13 


Problem  14 


Problem  17 


Problem  18 
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Abstract 


In  a  military  theater,  a  threat  region  can  be  modelled  as  a  convex  set  in 
B?.  A  mission  plan  for  a  cruise  missile  is  a  path  from  its  launch  location 
to  its  target  location  which  avoids  the  threat  regions,  if  possible,  and  does 
not  violate  a  turn  angle  restriction.  A  strike  plan  for  a  cruise  missile  strike 
involving  n  missiles  is  a  set  of  n  unique  mission  plans,  all  of  which  avoid 
the  threat  regions,  if  possible.  This  manuscript  presents  a  new  algorithm 
for  automatically  creating  a  strike  plan  in  the  presence  of  threats.  The  new 
algorithm  is  based  on  the  two-dimensional  geometry  of  the  problem  and 
mimics  a  strategy  which  a  human  operator  could  easily  apply.  Twenty  test 
cases  have  been  solved  and  graphical  displays  are  presented  for  each  test  case. 

Acknowledgement 

This  research  was  supported  in  part  by  the  Office  of  Naval  Research  under 


contract  number  N00014-96-1-0315. 


1  Introduction 


There  are  a  variety  of  missiles  which  are  classified  as  cruise  missiles.  In 
general,  they  carry  a  GPS  system,  a  terrain  following  system,  and  a  vision 
system  which  allows  them  to  use  TERCOM  maps  for  precise  navigation.  A 
mission  plan  for  a  cruise  missile  is  defined  as  a  path  composed  of  line  segments 
from  a  launch  site  to  a  target  site  which  does  not  violate  any  turn  restrictions. 
The  missile  is  programmed  with  the  waypoints  defining  the  mission  plan.  The 
GPS  navigation  system,  in  conjunction  with  a  TERCOM  map  system,  guides 
the  missile  to  its  intended  target  without  any  outside  intervention.  Once  a 
cruise  missile  is  launched,  its  mission  cannot  be  modified. 

Since  cruise  missiles  are  low  flying,  they  are  generally  not  susceptible  to 
attack  by  enemy  SAM  sites.  However,  they  are  vulnerable  to  anti-aircraft 
artillery  (AAA),  if  a  ground  crew  has  their  weapons  aimed  at  exactly  the 
location  over  which  the  missile  will  fly.  For  example,  if  two  missiles  fly  the 
same  mission  with  a  ten  minute  time  interval  between  flights,  the  ground 
crews  along  the  flight  path  have  an  excellent  chance  of  destroying  the  second 
missile.  The  current  strategy  used  to  defeat  AAA  is  to  vary  the  missions 
so  that  no  two  missiles  fly  identical  missions.  Varying  a  mission  by  a  few 


hundred  yards  to  either  the  right  or  left  is  usually  sufficient  to  defeat  AAA. 

A  mission  plan  is  for  a  single  missile  and  a  collection  of  mission  plans  is 
called  a  strike  plan.  During  recent  conflicts  in  the  Persian  Gulf,  cruise  missile 
strike  plans  involved  approximately  25  missiles.  These  strikes  generally  occur 
early  in  the  morning,  involve  a  two  to  three  hour  time  window,  and  require 

approximately  25  unique  mission  plans. 

In  a  previous  investigation  [6],  a  new  algorithm  was  developed  to  auto¬ 
matically  generate  a  cruise  missile  mission  plan.  The  plan  (path)  is  composed 
of  line  segments  from  a  given  origin  (missile  location)  to  a  given  destination 
(target  location)  which  avoids  threat  regions  modelled  as  circles  and  convex 
polygons.  The  algorithm  is  based  solely  on  the  geometry  of  the  problem  and 
attempts  to  mimic  the  trial-and-error  strategy  a  human  mission  planner  could 
apply  to  this  graphical  problem.  The  basic  strategy  for  creating  straight  line 
paths  around  obstacles  is  to  move  along  the  edges  of  triangles  which  enclose 
the  obstacles.  A  heuristic  strategy  is  applied  to  determine  the  stopping  point 
along  the  two  sides  of  the  enclosing  triangle  and  branch- and-bound  is  used 
to  ensure  short  paths.  Unlike  optimization  problems,  the  mission  planning 
problem  as  described  in  [6]  has  no  mathematical  objective  function.  The 
objective  is  to  find  short  feasible  missions  with  only  a  few  segments.  The 


objective  of  this  study  is  to  present  an  algorithm  to  automatically  create  a 
strike  plan  involving  25  unicjue  missile  mission  plans. 

1.1  Description  of  the  Problem 

A  mission  plan  for  a  single  missile  is  a  set  of  connected  line  segments  from 
an  origin  to  a  destination  which  does  not  violate  a  turn  angle  restriction  and 
does  not  intersect  any  threat  region,  if  possible.  The  threat  regions  are  mod¬ 
elled  as  circles  and  convex  polygons.  A  strike  plan  is  a  collection  of  feasible 
unique  mission  plans.  A  feasible  strike  plan  involving  two  missiles  located 
at  M  with  target  T  is  illustrated  in  Figure  1.  One  mission  is  defined  by  the 
line  segments  {MA,AT)  with  turn  angle  /.CAT  and  the  other  is  {MB,BT) 
with  turn  angle  ITBD.  Another  strike  plan  with  the  same  two  missiles  and 
targets  is  illustrated  in  Figure  2.  Note  that  the  two  missions  are  similar  but 
unique.  Missions  which  differ  by  a  few  hundred  yards  to  the  left  or  right  are 
generally  different  enough  to  defeat  AAA.  A  feasible  strike  plan  having  three 


missiles  is  illustrated  in  Figure  3. 


Figures  1,  2,  3  about  here 


Using  the  ideas  and  definitions  presented  in  [6],  the  problem  of  finding  a 
feasible  mission  plan  may  be  stated  as  follows: 

Given  two  points  M  and  T  e  turn  angle  restriction  Z, 

and  K  convex  sets  i2i, . . . ,  Rki  find  a  set  of  line  segments 
[Yi, F2],  [>2, 5^3],  •  •  • ,  [>s,  Ys+i]  with  Fi  =  M  and  ys+i  =  T 

C 

such  that  no  line  segment  intersects  the  interior  of  any  Ri,..  .,Rk 
and  no  turn  angle  exceeds  Z. 

A  strike  planning  problem  involving  L  missiles  is  to  develop  L  unique  feasible 
mission  plans.  For  a  problem  having  L  missiles,  there  are  an  infinite  number 
of  feasible  strike  plans,  and  our  algorithm  attempts  to  create  a  plan  having 
three  essential  characteristics:  the  distance  for  each  of  the  mission  plans  is 
short,  each  mission  uses  a  relatively  small  number  of  waypoints,  and  no  turn 
angle  restriction  is  violated. 


1.2  Survey  of  the  Literature 

Several  investigations  for  mission  planning  cover  the  area  with  a  grid  graph 
and  apply  a  shortest  path  procedure  to  find  paths  through  the  grid  (see 
[1,  5, 16, 17,  27,  28]).  This  strategy  converts  a  problem  that  is  fundamentally 
continuous  into  one  that  is  discrete.  The  resulting  paths  must  be  adjusted 
so  that  there  are  only  a  few  line  segments,  and  it  is  difficult  to  include 
the  turn  angle  constraint  using  this  strategy.  We  believe  that  the  strategy 
implemented  in  [6],  which  exploits  the  geometric  properties  of  the  problem, 
is  a  superior” approach.  We  know  of  no  investigations  which  involve  multiple 
missiles  in  the  same  region. 

The  robotics  literature  on  path  planning  in  the  presence  of  obstacles  is 
also  related  to  the  problem  of  planning  missions  for  cruise  missiles.  Algo¬ 
rithms  for  generating  shortest  paths  for  mobile  robots  are  generally  based  on 
finding  the  shortest  path  in  the  tangent  graph  (or  visibility  graph)  consisting 
of  the  tangents  between  obstacles  (see  [7,  10,  11,  12,  13]).  When  applied 
to  mission  planning,  this  approach  can  result  in  numerous  very  short  seg¬ 
ments  and  these  ideas  do  not  easily  adapt  to  curved  surfaces.  We  know  of 
no  investigations  involving  path  planning  for  multiple  mobile  robots. 


2  The  Strike  Planning  Algorithm 


The  K  convex  threat  regions  are  denoted  by  i?i, . . . ,  Rk,  each  with  respective 
center  Xi,...,  Xk-  If  threat  region  Ri  is  a  circular  region,  it  has  an  associated 
radius  r,-.  If  threat  region  Ri  is  defined  as  a  convex  hull,  it  will  have  an 
associated  set  of  gi  generator  points  {G’j, . . . ,  G^.},  in  which  case  the  center 
Xi  is  defined  as  the  average  of  the  generator  points  so  that 

The  L  missiles  involved  in  the  strike  are  assumed  to  be  spread  out  among 
P  origin-destination  pairs  of  missile  locations  denoted  by  and 

target  locations  denoted  by  ti,  . . . ,  Tp,  with  Ai, . . . ,  Ap  missiles  assigned  to 
the  respective  origin-destination  pairs.  Thus  A,-  missiles  will  be  fired  from 
missile  location  fii  at  target  location  r,-.  Furthermore  L  =  Ai  -| - -1-  Ap. 

If  a  clear  path  (not  intersecting  any  threat  region)  is  available  from  missile 
location  m  to  target  location  r,-,  the  strike  planning  algorithm  generates  A,- 
two-segment  paths  {YiY^Ys, . . .  ,YiY2'Y3},  where  Yi  =  /i,-,  Y3  =  r,-,  and 
the  points  Y^,---,  Y2'  are  spread  out  equally  spaced  about  the  midpoint 
^{Yi  -}-  is)  along  the  line  perpendicular  to  line  seqment  Y1Y3  and  passing 


through  that  midpoint.  The  topmost  and  bottommost  sets  of  paths  in  Figure 
7  were  generated  in  this  manner. 

If  a  clear  path  is  not  available  from  missile  location  fii  to  target  location 
T,-,  the  strike  plannning  algorithm  makes  use  of  the  circumscribed  triangle 
algorithm  to  generate  a  set  of  A,-  paths  from  missile  location  m  to  target 
location  Ti,  which  will  have  good  separation. 

The  circumscribed  triangle  algorithm  starts  at  M  (a  missile  location)  and 
constructs  line  segments  by  moving  from  a  given  segment  point  tangentially 
to  the  circumference  of  nearby  threat  regions,  until  a  suitable  straight-line 
path  from  a  segment  point  to  T  (a  target  location)  is  found.  This  basic 
strategy  is  embedded  in  a  branch-and-bound  framework,  so  that  a  minimal 
distance  path  of  this  restricted  type  not  violating  turn  angle  limits  will  be 
obtained. 

When  the  strike  planning  algorithm  makes  use  of  the  circumscribed  tri¬ 
angle  algorithm  with  a  given  origin-destination  pair  (/ii,T,),  it  generates  a 
sequence  of  A,-  single  missile  paths  with  M  =  //,•  and  T  =  Ti.  After  each  sin¬ 
gle  missile  path  is  generated,  all  threat  regions  which  are  tangent  to  any  of 
the  path  segments  are  expanded  slightly.  For  a  circular  region,  the  radius  is 


slightly  increased.  For  a  convex  hull  region,  all  generator  points  Gj, . . . ,  G],. 
are  projected  radially  slightly  away  from  the  center  Xi.  This  expansion  in¬ 
sures  that  on  the  next  single  missile  path  generation,  the  same  path  will 
not  be  followed.  In  many  Ccises,  the  next  path  generated  will  have  the  same 
general  shape  but  be  deflected  slightly  away  from  the  previous  path  at  the 
points  of  tangency  to  the  threat  regions  encountered.  An  example  of  such  a 
set  of  paths  is  the  topmost  set  in  Figure  4.  In  some  cases,  one  or  more  paths 
may  take  a  slightly  different  route.  An  example  of  such  a  set  of  paths  is  the 
topmost  set  in  Figure  23.  This  indicates  that  other  good  paths  are  close  in 
distance  to  the  minimal  length  restricted  path  found.  After  all  A,-  paths  have 
been  generated  in  this  manner  for  a  given  origin-destination  pair  all 

modified  threat  regions  are  returned  to  their  originally  specified  shapes. 


3  Empirical  Analysis 

Using  a  grid  size  of  780  x  780,  five  groups  of  test  problems  were  randomly 
generated  and  strike  plans  were  developed.  The  four  problems  in  each  group 
have  the  same  threat  regions  composed  of  a  combination  of  convex  polygons 
aaid  circles,  but  the  origins  ajid  destinations  of  the  25  missiles  may  be  dif- 


ferent.  All  missiles  are  located  on  the  left  boundary  with  all  targets  located 
on  the  right  boundary.  Each  of  the  twenty  test  problems  were  solved  with 
the  strike  planning  algorithm  and  the  solutions  are  displayed  in  Figures  4 
through  23.  The  plots  were  obtained  using  Tcl  and  Tk  Solver  [23],  which  is 
part  of  our  experimental  computational  package.  In  all  cases,  the  algorithm 
produced  a  strike  plan  having  the  desired  characteristics. 

4  Summary  and  Conclusions 

Using  the  mission  planning  algorithm  described  in  [6]  as  the  fundamental 
computational  engine,  this  manuscript  presents  a  new  algorithm  to  auto¬ 
matically  generate  a  strike  plan  for  a  cruise  missile  strike.  A  strike  having 
25  missiles  is  composed  of  25  unique  mission  plans,  each  of  which  consists 
of  a  path  from  its  location  to  its  target.  Even  though  two  missiles  located 
on  the  same  ship  or  submarine  have  the  same  target,  their  missions  must  be 
different  enough  so  that  they  are  not  vulnerable  to  AAA. 

Like  the  mission  planning  algorithms  presented  in  [6],  this  algorithm  ex¬ 
ploits  the  geometry  of  the  problem  description  and  mimics  the  trial-and-error 
approach  currently  in  use  by  the  Navy’s  mission  planners.  Since  there  is  no 


single  criteria  for  evaluating  a  strike  plan,  the  merit  of  the  solutions  can 
only  be  judged  qualitatively  by  examination  of  the  graphical  display  of  the 
solution.  We  believe  that  the  twenty  displays  presented  in  this  manuscript 
provide  a  convincing  case  that  our  strike  planning  algorithm  provides  good 
strike  plans. 
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